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The large mass of the top quark, compared to that of 
the other fundamental particles, gives it a special position 
within the standard model (SM). Since its discovery [l|, 
the top quark has played an important role in theoret- 
ical extensions beyond the standard model (BSM) 0]. 
Recently, renewed interest has been directed toward 
searches including top quark final states for BSM physics 
due to discrepancies reported in the ti forward-backward 
asymmetry [g, 0]- Moreover, the most recent search 
for resonant ti production from DO Q reports an ap- 
proximately 2(7 excess of events at resonant-mass values 
around 950 GeV/c 2 . Many BSM theories predict 
heavy resonances that add a resonant component to the 
SM ti production mechanism. 

Top quarks decay via the weak interaction, nearly al- 
ways into a W boson and a b quark. The W bosons 
then decay into lighter fermion-antifermion pairs (ll| . A 
leptonic decay into a charged lepton and a neutrino oc- 
curs 32% of the time while a hadronic decay into an up- 
type quark and a down-type quark occurs the remain- 
ing 68% of the time. We search for resonant produc- 
tion of top quark pairs followed by decays into a final 
state with one lepton and multiple jets, where one of the 
W bosons decays leptonically (to either an electron or 
a muon plus a neutrino) and the other W boson decays 
hadronically. This semi-leptonic channel features a dis- 
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tinctive final state due to the presence of a charged lepton 
and has a branching ratio of 29%. 

we do not 



Unlike previous searches at CDF [12H15 



apply constraints based on the presence of top quarks 
in the event. While we focus the discussion on ti reso- 
nances, we construct the top-antitop mass M t i used as 
a final search discriminant by taking the invariant mass 
of all objects (lepton, jets, and missing Et) in the event 
including those that may not originate from top quark 
production. Other than the event selection defined be- 
low, which provides a sample primarily composed of ti 
events, there are no requirements that the event be con- 
sistent with ti production. This results in a more general 
search that is sensitive not only to ti, but also to any 
heavy narrow resonance decaying into a final state with 
a W boson and three or more jets with one or two jets 
originating from a b quark. 

As a benchmark model, we consider a specific SM ex- 
tension, topcolor-assisted technicolor [l6j]. This model 
explains the large mass of the top quark through the in- 
troduction of new strong dynamics and also predicts a 
vector particle (Z' boson), which couples primarily to 
the third generation of quarks and has no significant 
couplings to leptons. The existence of a narrow- width 
Z' boson resonance (IV = 0.012 Mz>) decaying to ti 
pairs, using the leptophobic topcolor model [l7j. has been 
searched for both by the CDF and DO @, EE [1^1 

experiments at the Tevatron, and also by the ATLAS 
and CMS 
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|22| experiments at the LHC. For reso- 
nance searches at the highest masses, the LHC experi- 
ments have superior sensitivity to the Tevatron due to 
the higher center-of-mass energy. However, in the lower 
mass regions (mz< < 750 GeV/c 2 ) the Tevatron experi- 
ments have competitive sensitivity in searches for parti- 
cles produced in gg-initiated states, such as the Z' boson. 
While the production rate for the main background from 
SM ti production is approximately 25 times larger, no 
valence antiquarks are available in the LHC pp collisions, 
so the signal production rate increases by a smaller fac- 
tor relative to the pp collisions of the Tevatron (between 
four and eight depending on the signal mass hypothesis). 

The collision events discussed in this Letter were pro- 
duced at the Tevatron pp collider at a center-of-mass en- 
ergy of 1.96 TeV and were recorded by the CDF II de- 
tector [23| . The data sample corresponds to the full data 
set of the Tevatron, which comes from an integrated lu- 
minosity of 9.45 fb _1 . The CDF II detector consists of 
high-precision tracking systems for vertex and charged- 
particle track reconstruction, surrounded by electromag- 
netic and hadronic calorimeters for energy measurement, 
and muon subsystems outside the calorimeter for muon 
detection. CDF II uses a cylindrical coordinate system 
with azimuthal angle <p, polar angle 6 measured with re- 
spect to the positive z direction along the proton beam, 
and the distance r measured from the beamlinc. The 
pseudorapidity, transverse energy, and transverse mo- 



mentum are defined as ?] = — In [tan(|)l , Et — EsinO, 
and pt = psin#, respectively, where E and p are the en- 
ergy and momentum of an outgoing particle. The miss- 
ing transverse energy St is defined by St — — Yli Er^ii 
where hi is a unit vector perpendicular to the beam axis 
that points to the iih calorimeter tower (St = \St | ) • 

The event selection and background estimation meth- 
ods summarized below closely follow those that were em- 
ployed in the observation of single top quark produc- 
tion [24[ and in the search for the Higgs boson in the 
WH — > Ivbb final state [25[. The main difference with 
respect to the current search is the jet multiplicity re- 
quirement. 

The data were collected using online event selections 
(triggers) requiring one of the following energetic-lepton 
signatures: a high transverse momentum (pt) electron 
candidate, a high-px muon candidate, or large St- Sig- 
nificant St can be produced when the neutrino from a 
leptonic W boson decay escapes detection. 

Candidate events are selected by requiring a lepton 
candidate with p e T > 20 GeV/c, St > 20 GeV, and three 
or more jets with | ^7 1 < 2.0 and Et > 20 GeV after cor- 
recting the jet energies for instrumental effects 2(| 27 1. 
One or two jets must be identified as being likely to have 
originated from a b quark according to the SECVTX [28| 
algorithm. This algorithm searches in the jet for a sec- 
ondary vertex which results from the displaced decay of 
B hadrons. Events are rejected if more than one iden- 
tified lepton is reconstructed, or if they are kinemat- 
ically inconsistent with leptonic W boson decays [29| . 
Events with severely misreconstructed jets or leptons are 
removed based on angular correlations between the jet or 
lepton candidate and the Jr- 

Models for background processes are derived from a 
mixture of simulation and data-driven techniques (24J. 
Important backgrounds in this final state include SM ti 
production and other processes that include a W or Z 
boson in association with three or more jets. The events 
can include true 6-quark jets, as in W boson + bbj events, 
or jets that have been misidentified as 6-quark jets, such 
as in W boson +ccj and W boson +jjj events, where 
j refers to jets not originating from heavy-flavor quarks. 
Multijet events without W bosons also contribute to the 
sample composition. Additional small background con- 
tributions come from Z boson production with additional 
jets, diboson production, and single top quark produc- 
tion. 

The expected rate for the SM ti background is taken 
to be 7.04 ± 0.50 pb as calculated at ncxt-to-next- 
to-leading order using MSTW 2008 parton distribution 
functions [3l| . In order to predict the acceptance for non- 
resonant SM ti events and their kinematic distributions, 
we use a sample of Monte Carlo (MC) events generated 
using POWHEG 1321 and assuming a top quark mass of 
172.5 GeV/c 2 [331 with parton showering provided by 
pythia v6.2 [34| followed by simulation of the CDF II 
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Process 


3- jet events 


> 4-jet events 


tt 


1925 ± 204 


2565 ± 271 


W/Z boson + jets 


2281 ± 607 


569 ± 189 


Multijets 


147 ± 60 


126 ± 104 


Total background 


4354 ± 872 


3260 ± 563 


Observed 


4254 


3049 



TABLE I: Summary of the background prediction and ob- 
served data for three-jet and four- or- more- jet events. The 
uncertainties include statistical and systematic contributions. 



the neutrino is estimated from the Mt- However, be- 
cause the z-component of the momenta of the scattering 
partons from the pp collision is unknown, the final-state 
reconstructed energy need not be balanced in the z di- 
rection. The longitudinal component of the neutrino mo- 
mentum (p z ) is determined by solving = (p l +p") 2 , 
where Mw, p l , and p v are the W boson mass, the lepton 
momentum, and the neutrino momentum, respectively. 
The smaller solution of the resulting quadratic equation 
is chosen for the p z of the neutrino. If there is no real 
solution we set neutrino p z = 0. This approach is found 
to select the correct p z of the neutrino in about 70% of 
simulated tt events. 



detector [35|, |36[ . The detection efficiency predicted by 
the MC is corrected based on measurements using data 
for lepton identification, trigger efficiencies, and 6-jet tag- 
ging efficiencies. The normalization for the W boson + 
jet processes is obtained from a fit to the $t distribu- 
tion before the 6-tagging requirement is applied. The 
background from events with mistakenly 6-tagged light- 
flavor jets, W boson for example, is estimated by 
measuring the rate of such mistags in multijet data [28 1 
and applying this rate to the W boson + jets data sam- 
ples before b tagging. The contribution from true heavy- 
flavor production in the W boson + jets event sample is 
determined from measurements of the heavy-flavor event 
fraction in a W boson + 1 jet sample that is independent 
of the sample used in the resonance search. We model 
the kinematic distributions of W boson + jets events us- 
ing a combination of ALPGEN [37| matrix-element gen- 
eration and PYTHIA parton showering. The rate of the 
QCD multijet background is obtained from a fit to the 
$t distribution, where the QCD multijet background is 
modeled using a sample of collision events in which one 
of the lepton identification requirements is inverted to 
obtain an enriched sample of QCD multijet events. 

The background predictions are summarized in Tablc[TJ 
In this table and the following figures we have divided the 
sample into events that include three jets and events that 
include four or more jets. For the statistical interpreta- 
tion of the data we further subdivide the events based on 
the number of 6-tagged jets (one or two b tags) and based 
on the lepton type (lepton types that can be directly 
identified by the trigger, or leptons in events selected 
with the -E^-based trigger), yielding eight independent 
channels used to search for a resonance in the M t j distri- 
butions. The sensitivity of the search benefits from this 
subdivision because the search subchannels have differ- 
ent background compositions, signal-to-background ra- 
tios, and invariant mass resolutions. 

We use the invariant mass of all reconstructed objects 
in the event to discriminate between SM background and 
Z' boson signal events. For each event we calculate M t j 
using the momenta of the three or more jets, the charged 
lepton, and the neutrino. The transverse momentum of 









Acceptance 




o z ,BR(Z' -> tt) 


[GeV/c 2 ] 


3 jets 


[%] 


>4 jets [ 


%] 


[Pb] 


350 


2.75 


± 


0.15 


3.68 


± 


21 


8.91 


400 


2.82 


± 


0.16 


3.83 


± 


21 


12.3 


450 


2.35 


± 


0.13 


3.47 


± 


19 


8.24 


500 


2.29 


± 


0.13 


3.59 


± 


20 


5.53 


550 


2.16 


± 


0.12 


3.63 


± 


21 


3.51 


600 


1.93 


± 


0.11 


3.58 


± 


21 


2.30 


650 


1.71 


± 


0.10 


3.48 


± 


20 


1.43 


700 


1.52 


± 


0.09 


3.24 


± 


19 


0.917 


750 


1.37 


± 


0.09 


2.99 


± 


18 


0.566 


800 


1.19 


± 


0.08 


2.75 


± 


17 


0.355 


850 


1.10 


± 


0.07 


2.46 


± 


15 


0.208 


900 


0.96 


± 


0.06 


2.24 


± 


14 


0.134 


950 


0.92 


± 


0.06 


2.03 


± 


12 


0.080 


1000 


0.87 


± 


0.06 


1.91 


± 


12 


0.049 


1100 


0.86 


± 


0.06 


1.84 


± 


11 


0.017 


1200 


1.05 


± 


0.07 


2.01 


± 


12 


0.006 



TABLE II: Acceptances and predicted cross sections times the 
branching ratio for the leptophobic topcolor Z' boson signal 
hypotheses. 



For the benchmark model, the Z 1 boson cross sections 
times branching fraction are based on leading-order pre- 
dictions from Ref. [38| with an additional scaling factor 
of 1.3 applied to account for next-to-leading-order (NLO) 
effects [39j . Signal Z' boson events are modeled with sim- 
ulated events generated by PYTHIA in order to study the 
signal acceptance and to predict the M t i distributions. 
Table [TTJ shows the selection efficiencies and cross sec- 
tions for Z' boson events after the final event selection 
for each mass hypothesis considered in the analysis. 

A total of 4254 (3049) events survive the selection cri- 
teria for the three-jet (four-or-more-jet) category. The 
SM tt contribution is estimated to be 43% (78%) for 
three-jet (four-or-more-jet) events. The remaining events 
are contributed primarily from the W boson + jet and 
QCD multijet processes plus a potential signal contribu- 
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FIG. 1: Reconstructed M t f for the three-jet events (left) and four-or-more-jet events (right). The full distribution is on a linear 
scale while the high-mass tail in the inset is drawn on a logarithmic scale. The background expectation is normalized to the 
best fit from the data. The red histogram shows the expectation for a 600 GeV/c 2 leptophobic topcolor resonance normalized 
to the predicted cross section. 



tion from Z' boson events. The M t j distributions for the 
background model and events observed in the data are 
shown in Fig. [TJ The M t i distribution for the Z 1 bo- 
son signal for the 600 GeV/c 2 mass hypothesis is also 
included in Fig. [TJ 

We calculate a Bayesian credibility level (C.L.) limit 
on resonant tt production for each mass hypothesis based 
on the binned observed M t i spectrum using the combined 
likelihood which includes the priors, tt(6, on the system- 
atic uncertainties, 9: 



function, multiplied by the 9k priors, jr(9k), is then inte 



grated over 9k including correlations [111 ]: 



,-9*/2 



C(R, s, b\n, 9) x n(9) = jj fl C ~ZTT X II 
i=l j=l ni 3' fc=l 

(1) 

In this expression, the first product is over the number of 
channels Nq, and the second product is over histogram 
bins containing rnj events. The predictions for the bin 
contents are /Lty = R x Sij(6) + bij{&) for channel i and 
histogram bin j, where represents the potential res- 
onant signal, bij is the expected background in the bin, 
and R is a scaling factor applied to the signal. 

Systematic uncertainties are parametrized by the de- 
pendence of Sij and b^ on 9. Each of the n sys components 
of 9, 9k, corresponds to a single independent source of 
systematic uncertainty. We account for correlations by 
allowing each parameter to have an impact on several 
sources of signal and background in different channels. 
Gaussian priors are assumed for the 9k, which are trun- 
cated so that no prediction is negative. The likelihood 



C'(R) = J C(R,s,b\n,9)n(9)d9. 



(2) 



We assume a uniform prior in R to obtain its poste- 
rior distribution. The observed 95% C.L. upper limit 

on R, R$? s , satisfies 0.95 = f£™ £'(n)dH. The ex- 
pected distribution of i?gs is computed in an ensemble 
of pseudoexperiments generated without signal. In each 
pscudoexperiment, values of the nuisance parameters are 
drawn from their priors. The median expected value of 
i?95 in this ensemble is quoted as the expected limit. This 
statistical procedure is repeated for each resonance-mass 
hypothesis from 350 GeV/c 2 to 1200 GeV/c 2 . 

We consider uncertainties that affect the normaliza- 
tion as well as uncertainties that affect the M t j distribu- 
tions. The same set of uncertainties on the dominant 
background (SM tt production) and the resonant sig- 
nal are considered: they arise from the uncertainty in 
the jet energy scale (JES) \21j, the ^-tagging efficiency, 
the luminosity measurement [40j, the lepton identifica- 
tion and trigger efficiency (2-6%), and the rate of initial- 
and final-state (IFSR) radiation from the parton shower 
model. The JES, £>-tag, and IFSR variations also affect 
the shape of the M t j distributions. The rate of produc- 
tion for events with a W boson and heavy-flavor jet (b or 
c) is assigned an uncertainty of 30% due to limitations in 
the calibration of the fraction of heavy-flavor jets in the 
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sample. Uncertainties on the renormalization and factor- 
ization scale used in the ALPGEN sample affect the shape 
of the M t f distributions from W boson + jets. The QCD 
multijet background normalization is assigned a 40% un- 
certainty due to statistical limitations from the fitting 
procedure and the definition of the multijet model Q . 

The resulting 95% C.L. upper limits on a(pp — > 
Z')BR(Z' — > tt) as a function of M t j are shown in 
Fig. [2] and Table Mil together with expected limits de- 
rived from pseudoexperiments that include the SM back- 
ground hypothesis only. A benchmark leptophobic top- 
color model is excluded at 95% C.L. for Z' boson masses 
smaller than 915 GeV/c 2 assuming the width of the res- 
onance is Y z 1 = 0.012 Mz'. In addition, the limits re- 
ported here can be applied to any resonance producing 
the same final state as long as the decay width is signif- 
icantly smaller than the reconstruction mass resolution 
(Tz' <C 0.15 Mz>), and the difference in the acceptance 
with respect to the values quoted in TablelUis taken into 
account. 




FIG. 2: Expected and observed upper limits on the produc- 
tion cross section times the branching ratio for the leptophobic 
topcolor Z boson in 9.45 fb _1 of CDF data. The dashed line 
is the median expected upper limit with the assumption of no 
signal, the black points are the observed limit, and the blue 
line is the cross section prediction for leptophobic topcolor Z' 
boson production. 

In conclusion, we have performed a search for a heavy 
resonance decaying into tt using the semi-leptonic decay 
channel in data from 9.45 fb" 1 of integrated luminosity. 
The data are found to be consistent with the background 
expectation and upper limits are set on the produc- 
tion cross section times branching ratio at the 95% C.L. 



M Z ' [GeV/c 2 ] 


Expected [pb] 


Observed [pb] 


350 


0.772 


0.687 


400 


0.575 


0.652 


450 


0.670 


0.585 


500 


0.520 


0.427 


550 


0.354 


0.530 


600 


0.245 


0.472 


650 


0.199 


0.269 


700 


0.159 


0.145 


750 


0.137 


0.112 


800 


0.115 


0.099 


850 


0.106 


0.103 


900 


0.097 


0.116 


950 


0.091 


0.118 


1000 


0.092 


0.129 


1100 


0.098 


0.132 


1200 


0.134 


0.166 



TABLE III: Expected and observed 95% C.L. upper limits on 
the production cross section times the branching ratio for a 
narrow tt resonance, given as a function of Z' boson mass. 



For a specific benchmark model (leptophobic topcolor), 
we exclude Z 1 bosons with masses up to 915 GeV/c 2 . 
For masses smaller than approximately 750 GeV/c 2 , this 
search yields the most constraining limits to date on qq- 
produccd narrow tt resonant states in the semi-leptonic 
decay mode. 
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